Magneto shot noise in noncollinear diffusive spin-valves 
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We develop a semiclassical Boltzmann-Langevin theory of the spin polarized shot noise in a 
diffusive normal metal spin-valve connected by tunnel contacts to ferromagnetic reservoirs with 
noncollinear magnetizations. We obtain basic equations for correlations of the fluctuating spin- 
charge distribution and current density matrices by taking into account the spin-flip processes and 
precession of the spin accumulation vector in the normal metal. Applying the developed theory to a 
two terminal FNF structure, we find that for a small spin-flip strength and a substantial polarization 
of the terminals the shot noise has a nonmonotonic variation with the angle between magnetization 
vectors. While the shot noise is almost unchanged from the normal structure value for parallel 
configuration and increases well above the normal value for antiparallel configuration, it suppresses 
substantially at an intermediate angle depending on the ratio of the conductances of the N metal and 
the tunnel contacts. We also demonstrated pronounced effects of the polarization and the spin-flip 
scattering on the shot noise which reveals the interplay between relaxation and precession of the 
spin accumulation vector in the N metal. 

PACS numbers: 72.25.Ba, 72.25.Rb, 74.40.-|-k, 75.75.-|-a 



I. INTRODUCTION 

Spin-polarized transport in magnetoelectronic struc- 
tures has recently attracted an intense interest largely 
due to their important technological applications such 
as non-volatile magnetic random access memories 
(MRAMs), read heads for mass data storage and high 
sensitive magnetic sensorsiiSi^. The main effect is the 
(giant) magnetoresistance^ in magnetic multilayers and 
spin- valves, i.e. the remarkable decrease in the resistance 
when the orientation of the magnetization vectors of the 
ferromagnetic regions change from antiparallel to paral- 
lel. A spin-valve consists of two ferromagnetic leads as 
the spin injector and detector, connected through a nor- 
mal metal spacer which serves for the spin accumulation. 
Recent interest in transport with noncollinear (neither 
parallel nor antiparallel) magnetizations was stimulated 
by the spin-transfer magnetization torquc^'^, which is es- 
sential for the new devices such as the spin-flip2iS and 
spin torque transistorsiS. 

The most appropriate theoretical formalism for non- 
collinear ferromagnet-normal metal structures is the 
semiclassical circuit theory^ which is based on dividing 
of the structure into the reservoirs, nodes and the junc- 
tions and expressing the currents in terms of the scat- 
tering matrices of the junctions and the isotropic elec- 
tronic distributions inside the nodes and reservoirsii. For 
the noncollinear spin-transport the spin-charge current 
and the corresponding distributions inside the nodes and 
reservoirs, as well as the scattering matrices acquire the 
2x2 matrix characteristics of the Pauli spin space. The 
semiclassical Boltzmann equations provide the standard 
method to calculate the noncollinear distribution ma- 
trix in a continuous bulk medium in the diffusive limit. 
Magnetoelectronic circuit theory complemented by the 
Boltzmann diffusion equations have been widely used to 
study different magnetoelectronic DC effectai^. In spite 



of this there has been little attention toward the fluctu- 
ations of the noncollinear spin-polarized current. Time- 
dependent current fluctuations at low temperatures, the 
so called shot noise, provides valuable information about 
the transport process which are not extractable from 
the average currenti^ti^. In magnetoelectronic struc- 
tures, in which the spin of electrons plays an essential 
rule, the shot noise is expected to contain spin-resolved 
information, including spin-dependent correlations and 
spin accumulation and relaxation. This together with 
the importance of the noise in magnetoelectronic devices 
in view of applications motivate studying of the spin- 
polarized current fluctuations. 

The theory for current fluctuations in disorder con- 
ducting systems can be formulated by including the 
Langevin sources of fluctuations due to the random scat- 
tering from disorders and the fluctuations of the dis- 
tribution in the semiclassical Boltzmann equationsi^. 
The resulting Boltzmann-Langevin (BL) formalism has 
been used to calculate the shot noise in many meso- 
scopic structures and its results are widely believed to 
be identical with an ensemble average of the exact quan- 
tum results^?.. One important example is the universal 
one-third suppression^"'''^'' of the shot noise in a diffu- 
sive normal metal, which was correctly predicted within 
the BL approachiii^. Recently spin-polarized shot noise 
has been studied theoretically in normal metals con- 
nected by the ferromagnetic terminals with coUinear 
magnetization aJ^^'^"-^^'^'^ . In Ref. |23 a semiclassical BL 
theory of the collinear spin-polarized current fluctuations 
in a diffusive normal metal was developed. It was found 
that in a multi-terminal spin-valve structure the shot 
noise and the cross correlations measured between cur- 
rents of two different ferromagnetic terminals can devi- 
ate substantially from the unpolarized values, depend- 
ing on the spin accumulation and the spin-flip scatter- 
ing strength. On the other hand only a few works22i24 



have been devoted to shot noise in noncolhnear systems. 
Very recently Braun et.aP'^- showed that the frequency- 
dependent shot noise in a quantum dot spin-valve with 
noncolhnear magnetizations of ferromagnets in the pres- 
ence of an external magnetic field can be used to detect 
single-spin dynamics in the quantum dot. 

Tserkovnyak and Brataas^'' extended the circuit the- 
ory to obtain the Landauer-Biittiker (LB) formula for the 
shot noise of a FN contact with noncolhnear magnetiza- 
tion. They analyzed the angular dependence of the shot 
noise in a N metal node connected to two F reservoirs for 
different types of junctions between the reservoirs and 
the node. However, the effect of the spin-flip scattering 
inside the N metal as well as the spatial variation of the 
noncolhnear spin accumulation vector, which are essen- 
tial in diffusive spin-valves with wire geometry, are dis- 
regarded so far. These effects require an extension of the 
BL method to include the diffusion and the relaxation of 
the noncolhnear spin accumulation inside the diffusive N 
metal. To our knowledge there is no BL theory for shot 
noise of a noncolhnear spin-polarized transport. The aim 
of the present work is to develop such a theory. 

In this paper we develop a semiclassical BL equation 
for the noncolhnear spin-polarized current fluctuations in 
the presence of the spin-flip scattering. We obtain the ba- 
sic diffusion equations for the fluctuating spin-charge dis- 
tribution and current matrices which allows to calculate 
the mean current and the correlations of the correspond- 
ing fluctuations in a diffusive normal metal connected by 
tunnel contacts to several F reservoirs. The developed 
BL equations are supplemented by the the generalized 
LB formula for the shot noise at the contact points^. 

To illustrated the main behavior of the noncolhnear 
shot noise in magnetoelectronic systems we apply the 
developed BL equations to calculate the Fano factor, de- 
flned as the ratio of the noise power to average current, in 
a two terminal FNF structure. The noncolhnear orienta- 
tion of the magnetizations causes a precession of the spin 
accumulation vector through the N metal, which in pres- 
ence of the spin-flip scattering is associated with a damp- 
ing as a function of the distance from the F reservoirs. 
For a small spin-flip intensity and a flnite polarization of 
the tunnel contacts the precession of the spin accumula- 
tion results in a nonmonotonic angular dependence of the 
Fano factor. For a parallel conflguration Fano factor is 
almost the same as the the normal state value but for an- 
tiparallel conflguration it increases well above this value 
due to a large spin accumulation in the N metal. We flnd 
a substantial decrease of the Fano factor at an interme- 
diate orientation determined from the conductances of 
the N metal and the tunnel contacts. Introducing the 
spin-flip scattering as well as decreasing the spin polar- 
ization diminish the nonmonotonic behavior. We present 
a full analysis of the magneto shot noise which demon- 
strates the effects of the spin-flip induced relaxation and 
the precession imposed by noncolhnear magnetizations. 

The paper is organized as follows. Section II devotes 
to introduce the BL equation for noncolhnear transport 



and derive diffusion equations. In section III we calcu- 
late all possible correlations of intrinsic fluctuations and 
Langevin sources. Section IV devotes to calculate fluctu- 
ations and average of the current in double barrier FNF 
system. We show results for charge current shot noise in 
section V and flnally give some conclusions in section VI. 



II. BOLTZMANN-LANGEVIN EQUATIONS 
FOR NONCOLLINEAR TRANSPORT 

In this section we develop a semiclassical BL formal- 
ism for fluctuations of the noncolhnear spin-polarized 
current in a diffusive normal metal connected through 
tunnel contacts to a number of ferromagnetic reservoirs. 
When the magnetization vectors of the reservoirs have a 
noncolhnear orientation, the spin accumulation vector in 
the normal metal will also have a noncolhnear direction 
with respect to the quantization axis (z). The semiclas- 
sical electronic distribution is then determined by a 2 x 2 
matrix in the Pauli spin space of the form 

/ /TT(k,r,i) /THk,r,i) \ 
/(k,r,t)- . (2.1) 

\ /iT(k,r,t) /^i(k,r,t) J 

The fluctuating matrix distribution function /(k, r, t) = 

/(k, r) + (5/(k, r, i), depends on the momentum k, the 
coordinate r and time t. It is convenient to expand 
/ into the charge and the spin vector distributions in 
terms of the 2x2 unit matrix and the Pauli matrices 
{1, <T = {ax, ay^Uz)}- The matrices for two distributions 
whose spin vector distributions are pointed in opposite 
directions take the form 

/±(k,r,t)-/e(k,r,t)i±^./;(k,r,t). (2.2) 

Here /c = (/^^ +/^^)/2 is the charge or the spin indepen- 
dent part of the distribution matrix and /s the spin dis- 
tribution vector whose z-component fsz = (/^^ — /''^^)/2 
is spin polarization along the quantization axis and 
the other two components fsx, fsy describe the spin- 
polarization oriented perpendicular to the quantization 
axis. Note that the two distributions /^ are transformed 
to each other by making the replacement /s -^ ^fs- 

The BL equation for the noncolhnear distribution ma- 
trix 1)2.1(1 is written as 

at 

where /""p^''^) is the collision integral for normal spin- 
independent impurity (spin-flip) scattering and ^""P(^f) 
the corresponding Langevin source of the current fluctu- 
ations. The matrix collision integrals are expressed as 
the following 

i"»p[/+] = y dk'i^™p(k,k')[/+(k') - /+(k)], (2.4) 

rV+J-] = Jdk'w^\k,k')[r{k') - /+(k)], (2.5) 



where W™P(k, k') is the rate of the impurity scattering 
in which an electron scatters from the state with mo- 
mentum k into k' without changing its spin state. The 
spin-flip scattering rate VF''^(k, k') describes the transi- 
tion from the state with momentum k and spin state 
\s, ±) to k' and the flipped spin state |s, =f), where \s, =f) 
are the up and down spin eigen states in the quanti- 
zation axis paraUel to the local spin polarization vector 
/s = l/s|s- In writing the expressions (|2.4() and (|2.5() 
we assumed that H/™P('^f)(k,k') = H/™P('^f)(k',k), which 
follows from the detail balance principle. We also ignored 
dependence of the scattering rates on the spin state of 
electron. 

The Langevin sources are expressed in terms of the 
fluctuating part of the matrix distribution as follows 

|™P= f dk'W'"^P{k,k')[6f+{k') - 5f+{k)], (2.6) 



e 



:sf 



dk'W''{k,k')[Sf-{k') - df+{k)]. (2.7) 



For a diffusive normal metal we apply the standard dif- 
fusive approximation. Assuming that all quantities are 
sharply peaked around Fermi level, the momentum vec- 
tor k is expressed in terms of the energy e and the di- 
rection of the Fermi momentum n. Then the following 
relations are hold for the elastic scattering of electrons 
by the normal impurities and spin-flip disorders: 

p^imp(sf)(i^^j^/) = —5{e ~ e')w'"^'P^'^\n,n'), (2.8) 
No 

where A^o is the density of states in the Fermi level. In 
the diffusive regime the electronic distribution is weakly 
anisotropic in the momentum space, and the distribution 
matrix can be expanded up to the linear term in n: 



/(n,e,r,t) = /o(£,r,t) -I- n.fi(e, r, i). 



(2.9) 



Here the anisotropic part of the matrix distribution is 
determined by the vector fi whose components are 2x2 
matrices in the spin space. The matrix current density 
is expressed as 

J(r,t) = J(r) +(5J(r,f) = {eNoVp/G) f deh, (2.10) 

where vp is the Fermi velocity. The isotropic part /o 
determines the fluctuating electrochemical potential: 

(^(r, t) = (^(r) -t- <5(^(r, t) = (1/e) J defo- (2.11) 

Replacing expansion 12. 9|) in Eq. 1)2. 3|l we obtain the dif- 
fusion equations for the current density and the potential 
matrices, which read 



1 



Tr{ip) 



(2.12) 



V.J^-^r^-^i)+<..C (2.13) 



2Tsf 

J = —a\7(p + f). 



In these equations a = e^{No/2)D is the conductivity 
and (.si = VDtsi is the spin-flip scattering length, where 
the diffusion constant is given by D = VpT/3. The im- 
purity and the spin-flip relaxation times are deflned by 
the following relations 



n 



= / dn'w™P(n,n')(n-n'), (2.15) 



— = 2 / dn'w'\n, n')(n - n'), (2.16) 

r-si J 

(2.17) 



1 1 



Tsf 



The matrix Langevin current fluctuations in Eq. (|2.14|l 
has the form 

f) = evFiNo/2)T /"(I^P + f f)nrfnd£. (2.18) 



Also, 

<3--C(r,t) = (eA, 



o/2)|rfednff'-^^ij, (2. 



19) 



is the matrix Langevin divergence term due to the non- 
conserved nature of the spin-flip scattering. These two 
fluctuating Langevin terms have zero average values. 

Note that we disregarded the explicit dependence of / 
on time, given by the term d/dt of d/dt in the left hand 
side of the BL equation, since we are interested in the 
zero frequency noise power. To obtain these equations we 
focused on the realistic limit when Timp <C Tgf and ignored 
the effect of spin- flip scattering on the conductivity of the 
diffusive normal metal. 



III. CORRELATIONS OF CURRENT 
FLUCTUATIONS 



Correlations of the Langevin sources of current 
fluctuations 



To calculate correlations between components of the 
fluctuating Langevin terms r) and C we refer again to the 
BL equation H2.3|l . But, now we choose the spin quanti- 
zation axis to be parallel to the local mean spin distribu- 
tion vector fs- Denoting the projection matrices along s 
and two perpendicular unit vectors s±i{i = 1,2), by Wj 
{— & ■ s) and Vsj^i respectively, the matrix distribution 
function can be expanded as 



1 



1 



r = 7t(/5+ + /5-)l ± o(/-> - frs-)^-^ ± 



(2.14) 



JE('5/^-,+ -'^^-,-)^^-' (3.1) 



where the fluctuating transverse components of the spin 
polarization vector have zero mean values. 



We substitute the fluctuating parts of the matrices 
(|3.1|) into the Langevin source terms in Eq. H2.3|) and 
obtain the result 



where 



^ 



iinp(sf ) 



]_/^imp(sf) 



^hnp(sf)^^^ 



fimp(sf)-, j^ 



1 



,^imp 



p(sf) 



E(C[ 



p(sf) 

+ 



cr^'v. 



(3.2) 



where 



^^(sf) ^ / di^isj^;Z(^,k)-sj^l(k,k')], 



jimp{si). 



fimp(sf) 



i(-)c 
imp(sf) , 



and we have introduced the current in individual 
impurity (spin-flip) scattering as J]™"^, (k, k') = 

VF™P("f)(k,k')/s,a(k)[l-/s,a'(k')] for component of 
spin in direction of s. 

Now we apply the central assumption of the BL ap- 
proach and assume that all scattering events are inde- 
pendent elementary processes and thus the correlations 
of the associated currents fluctuations obey the Poisso- 
nian relation: 



(Sf-"'"'^''^ (ki, k2, r, i)<5J™P^''' (kg, k4, r', t')) 

= SijSaia3^a2aiS{ki - k3)5(k2 - k4) 

x<5(r - r')6{t - i')C^fi!(ki,k2,r,t), (3.3) 



and 

('5Ca,a.(ki,k2,r,t)W|;,„3„,(k3,k4,r',i')) -0 (3.4) 

Using these relations we can calculate the correlations be- 
tween matrix elements of the Langevin sources (^""p^''^)), 
from that we obtain all the possible correlations between 
components of the vector matrices r) and ^. The results 
read 



yf^'tr-' i-'W - 



1 



(Tyf '^(r, t)r,^^ (r', t')) = -J,™<5(r - r')5{t t')a 

V 

'(3.5) 
^7^^"^*"^'Ens-,.-.(r), (3.6) 



{vf{r,t}C{r',t'))^0, 



(3.7) 



n.-^..'(r) = / de/,,.(e,r)[l - /,,.'(£,r)]. (3., 



Here fg^^{e,r) is the isotropic part of the up and down 
spin components of the mean distribution function with 
respect to the quantization axis parallel to s. So, if we 
know the mean matrix distribution function we can eas- 
ily calculate these correlations. It is a remarkable result 
that, just the longitudinal fluctuations of the distribution 
function has non-vanishing contribution to correlations 
and hence in shot noise. But, the transverse fluctuations 
of distribution function, actually make no contribution 
in the shot noise. 



B. 



Boundary conditions and correlations of 
intrinsic fluctuations 



The diffusion equations ii>.ll^liJ.14|) and Eqs. (|3.5I3.7|) 
are the extension of the BL equations obtained in Ref. 
I21I for the collinear spin-polarized shot noise to the non- 
coUinear magnetizations. These equations are a complete 
set of equations, which have to be implemented by the 
appropriate boundary conditions at the contacts of the 
normal metal to the ferromagnetic terminals. For a FN 
junction with a noncoUinear magnetization the expres- 
sion for the matrix current reads^ 

eF = G^'^u^ (/f - /n) u^ + G^^u^ (/f - /n) u^ 

(3.9) 



-G^^uVnu^ - (GT^)*uVnu 



where /n (/f) is the average matrix distribution func- 
tion in the normal (ferromagnetic) side of the con- 
tact, GTT.ii = ^ [M-E™ki;4l'] are the diagonal 
components of the conductance matrix describing cur- 
rents of electrons with spin parallel (up) and antipar- 
allel (down) to the magnetization vector of F, where 
M is the number of transverse modes in the contact 
channel and rjj^ (^Im) and r^^ (^im) are the spin up 
and down reflection (transmission) coefflcients in the ba- 
sis where the spin-quantization axis is parallel to the 
magnetization of F. The mixing conductance G^-'- — 

^[M - Y,nm dmirirn)*] > describcs the transfer of the 
spin-current perpendicular to the magnetization of F. 
Here u^^^'> = I [^ + (")''" ' ™] are projection matrices in 
the spin space in which m is the unit vector in direction 
of magnetization vector. 

To write an expression for the fluctuations of the 
matrix current (|3.9|l we follow Beenakker and Biittiker 
assumption^^ and consider that the time-dependent fluc- 
tuations have two contributions. The first contribution 
SI{t) is the intrinsic fluctuations of the matrix current, 
which is caused by randomness of the electron scatter- 
ings from the junction. This is the only relevant term for 
the shot noise of a single junction connecting two reser- 
voirs which are held at equilibrium. The second con- 
tribution can exist when the distribution functions in N 



metal and/or the adjacent reservoirs are fluctuating. For 
a contact between the N metal with fluctuating part of 
the matrix distribution (5/n and a F reservoir held at 
equilibrium, the expression for the matrix current fluctu- 
ations has the form 

eA/^ =51- gTTmTj/j^(i.)uT _ G^^u^5h{Y)u^ 

-G^^u'^5h{Y)u^ ~ (G^^Yu^8f^(y)v}. (3.10) 



The correlations of the intrinsic fluctuations for an ar- 
bitrary type of a FN contact with a noncoUinear mag- 
netization was calculated by Tserkovnyak and Brataas^i 
using a generalized LB approach. Here we use their re- 
sult and write the correlations between different elements 
of the intrinsic matrix current fluctuations C'"''" ^ = 
Jdt{Si°''^{t)SI"'>^'{0)), as follow 



C"'^"'^' ^ del J2 ^ 



^SiS2S3S4 



S1S2 '- 

V Ja'B\ J aB' \ J a' B \ J aB' 



(3.11) 



In this equation G"^"^ = W 'Y^r. 

QS1S-2S3S4, _ S—\1\J \^ 

•^ — h l^" l^nmr 



nm\' mn 



pSi („S2 

nm \ mn 



)*] and 

m'n) J 



are elements of the conductance matrix G and the shot 
noise matrix S respectively. The general expressions of 
these matrices are considerably simplified for a tunnel 
barrier, for which an expansion in terms of the small 
transmission coefficients Sr'^^ = 1 — r^^, leads to the 
following resultsSi 



^ 


G 


rv^ 






2 




G 


,S' = 






2 



l + p 1 
1 I- p 

2 + 2P2 + P2 + P 2 
2+P 2 2 2-P 

2+P 2 2 2-P 

2 2-P 2-P 2- 2P 



(3.12) 
(3.13) 



Where we have the total conductance G = G^^ + G^^, 
the polarization P = (G^T - G^^)/G, ReG^^ = G/2 and 
IniG^^ = for the tunnel contact. Thus by knowing 
the average distribution function matrices inside the con- 
nected N and F metals, we can calculate all possible cor- 
relations between the intrinsic current fluctuations. 

Now we have all requirements to calculate the shot 
noise in a ferromagnetic multi-terminal diffusive struc- 
ture. The first step is to solve the diffusion equa- 
tions p. 1212. 14(1 implemented by the boundary condi- 
tions which are temporal current conservation laws in 
the contacts to obtain the mean distribution matrix in 
N metal and the charge and spin currents fiuctuations in 
the end points of N metal in terms of the Langevin cur- 
rent fiuctuations and the intrinsic current fiuctuations 
at these points. Now, knowing the correlations between 
these fluctuations from Eqs. H3.5I3.7|I and H3.11|l . enable 
us to calculate the correlations of the charge and spin 
currents fluctuations. In the next section we demon- 



N 



v,=o 




Vo=V 



FIG. 1: A schematic view of the system. A diffusive normal 
metal wire connected to two ferromagnetic reservoirs through 
the tunnel junctions. Ferromagnetic reservoirs are in thermal 
equilibrium and magnetizations of them can orient in arbi- 
trary directions. 



strate the above developed formalism by studying the 
noncoUinear shot noise in a two terminal spin-valve. 



IV. ANGULAR SHOT NOISE OF A DOUBLE 
BARRIER F-N-F SYSTEM 

We consider a two terminal diffusive spin- valve system 
as shown in Fig^ Two ferromagnetic reservoirs Fi and 
F2 with a noncoUinear orientation 6 of the magnetization 
vectors are connected through tunnel contacts to a dif- 
fusive normal metal wire (N) of length L. The reservoirs 
Fi^2 are held at and V potentials, respectively. We con- 
sider a symmetric structure for which both of the tunnel 
contacts have the same conductance and shot noise ma- 
trices given by Eas. 1(3. 1213. 13|l . 

The mean distribution function matrix inside the N 
metal is obtained by solving the diffusion equation H2.12|l . 



The general solution of this equation can be written as 
the energy integral of the following relation: 



hi^) 



h + {f2-h){a + b-) 



1 



csinh(Aa;/L) + dcosh{\x/L) (/s - /i), (4.1) 



where the matrix distribution function at the terminals 
Fi_2 are determined by the Fermi distribution functions 

via /i = /FD(e)i and /2 = f-pTi{e-eV)i, and A = i/4f is 
a dimensionless parameter which measures the intensity 
of the spin-flip scattering. The scalars a, h and compo- 
nents of the vectors c = (01,02,03) and d = (^1,(^2,^3) are 
coefficients which have to be determined by the boundary 
conditions. The boundary conditions are the spin-charge 
current conservation rules which in the absence of spin- 
flip process at the tunnel contact i (i = 1, 2), read 



/f = o. 



(4.2) 



Where I^ and if are expressions for the mean matrix 
current through the zth contact which arc calculated by 
relations H3.9() and the diffusion equation H2.14|l . respec- 
tively. The mean charge current is determined by the 
coefficient h only: 



Tr[i) 



-2G^b, 



(4.3) 



where the normal conductance of the N wire is given by 
Gn = aA/L, A is the cross section area of the wire. 

To obtain correlations of the Langevin sources in Eqs. 
(|3.5I3.7|) we need the matrix distribution function in the 
quantization axis parallel to the mean spin distribution 
vector. This is achieved by diagonalizing the distribution 
function matrix to get fs,±ix). The result is 

k±{^) = fi + {f2-fi)[{a + b^) 

± |csinh(Aa;/L) + crcosh(Aa;/L)|j , (4.4) 

where /i^2 are the corresponding Fermi distributions in 
the terminals. Now we calculate the fluctuations of the 
matrix current A/^ at the contacts {x — 0, L) inside N 
wire. Using the equation for 63 (|2.14l) and Eq. 12.13|l we 
obtain 



A/^(0,L) = -a (p ds 



V0eO(.)(.)^i 



+V(/.,o(L)(a;) (<5(^-^^^^i 



+ 61(0,1), (4.5) 



where integration is over the surface of the N metal and 



61(0, L)= I dr 



i(rr(77).V)0cO(L)(a:)i 



,7-^i)-V)0,o(L)(:^) 



(4.6) 



Here the potential functions are defined as 



(j)coix) = 1 - T' 



(Pcl{x) = 
(t>so{x) = 
4'sl{x) = 



sinh[A(l - x/L)] 

sinh(A) 
sinh(Aa;/L) 
sinh(A) 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



Note that as a result of the spin-flip scattering the spin 
current and hence, its fluctuations are not conserved 
through the wire. 

Now we can apply the current conservation law for the 
fluctuations of the matrix currents at two contacts: 



A/^ 



A/f = 



(4.11) 



Replacing the expressions of the current fluctuations 
from Eqs. |I4.5I4.10|) and Eq. (|TTn|) in Eq. l|CT|l we 
obtain a system of eight linear equations whose solutions 
give the fluctuations of the matrix chemical potential 
6ip{0,L) at the points x = 0,L in terms of the intrin- 
sic matrix current fluctuations i5/i,2 and the Langevin 
matrix current fluctuations 61(0, L). Replacing this re- 
sult into Eq. (|3.1()(l the matrix current fluctuations are 
expressed in terms of the fluctuating Langevin sources of 
currents and the intrinsic current matrices as 



AT 



a(3 _ 



a'P' 






+C2',f3,Sr P (0) + Dl?p,5I'^ '^ (L) , (4.12) 



where A° 



B"'f„,, C ,a, and D"'f„, are the coefficients 

dp ' a p dp 



.a/3 

a'P' 

P, A, and 



which depend on G/Gj^s, P, A, and 9. From this result 
we can calculate the correlations of the charge current 
S — 2 J dt{AIcAIc) and the corresponding Fano factor 
F = S/2eIc- The resulting expressions are too lengthy 
to be written down here. In the next section we analyze 
the magneto shot noise, the dependence of the F on the 
the relative angle 0, for different values of the involved 
parameters. 



V. RESULTS AND DISCUSSION 

Let us start our analysis of the shot noise with the 
limiting case of negligible spin-ffip scattering, i.e. when 
A -^ 0, and highly resistive tunnel contacts G/G^ -^ 
0. In this limit we retrieve the results of Tserkovnyak 
and Brataas for charge current shot noise^i, that is, a 
monotonic variation of -F as a function of the relative 
angle of the magnetization vectors 6: 
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(5.1) 



For a finite G/G-^ the angular dependence of Fano fac- 
tor deviates from the above simple monotonic behavior. 
This is illustrated in Fig. where we have plotted F{9) 
for different values of G/Gt^ in the limit A ^ and the 
polarization P — 0.99. By increasing G/Gj^, F decreases 
below the value given by Eq. (|5.1|l and finds a nonmono- 
tonic dependence on 0. The nonmonotonic behavior is 
more pronounced for a value of G /G^ which depends on 
the polarization. The Fano factor develops a minimum at 
a finite angle 0^0, which depends on P and G/G-^. In 
the limit of G/G-^ 3> 1, F returns back to the constant 
normal state value 1/3, as it is expected. This holds for 
not perfectly polarized terminals P ^ 1. We found that 
for P = 1, Fano factor retains its nonmonotonic angular 
dependence in the limit G/Gn ^1- In this case F has 
a minimum well below 1/3 at a finite 9, a sharp peak at 
9 = 7: where it reaches one and a smooth peak at 9 = 
with the value 1/3. The corresponding behavior in the 
conductance of the system was studied in Ref. [l^- 

Next we study the effect of spin-flip scattering in the N 
metal. In Fig. Owe show dependence of Fano factor on 
A for different angles, when G = 1 and P = 1. As it can 
be seen from Fig. |3 a strong spin-flip scattering destroys 
any spin polarization and thus suppresses the dependence 
of the Fano factor on 9. In this limit (A ^ 1) P reduces 
to the normal state value independent of 9 
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(5.2) 



At a finite A the shot noise of different ^s are separated. 
The most strong variations happens for A — > 0. Increas- 
ing A from zero, P(A) passes through a minimum or max- 
imum, depending on 9, before reaching the normal state 
value at large As. In the case of Fig. |3the minimum and 
maximum happen for 9 < 7r/2 and 9 > n/2, respectively. 
While the minimum and maximum values of P are lo- 
cated at A = for 9 = and 9 = tt, they can happen at 
finite As for intermediate angles. 

To show the effect of contacts polarization P, we have 
plotted in Fig. 0] the Fano factor vs. P in the limit 
A — > 0, for different 9s. For P — ^ the Fano factor takes 
the normal state value (|5.2|l , irrespective of the value of 
9. Increasing P the dependence of P on is appeared. 
The Fano factor increases (decreases) monotonically with 
P for 9 = TT (9 = 0) and reaches a maximum (minimum) 
at P = 1. However for an intermediate angles F{P) 
can have a nonmonotonic variation with a maximum at 
P < 1. This is in contrast to the limit of the normal 
metal node geometry (G/Gn -^ 0), where a monotonic 
dependence on P was predicted (Eq. ((OJ). 



VI. CONCLUSIONS 

In this paper we have presented a semiclassical the- 
ory of the spin-polarized current fiuctuations in a diffu- 
sive normal metal which is connected by tunnel contacts 
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FIG. 2: Fano factor versus the relative angle cos(S) for difler- 
ent values of G/Gn (Gn = 1) in the limit of A ^ and for 
P = 0.99. 



G=1,P=1 




FIG. 3: Fano factor versus the spin-flip scattering intensity 
A = L/lgf for different values of relative angle 6 for a system 
with G = 1 and P = 1. 



to ferromagnetic terminals with noncoUinear magnetiza- 
tions. Based on the Boltzmann-Langevin approach, we 
have developed diffusion equations which allow for calcu- 
lation of the charge-spin distribution and current density 
matrices and the correlations of the corresponding fluctu- 
ations in noncoUinear multiterminal systems. Our theory 
takes into account the precession of the noncoUinear spin 
accumulation as well as the spin-flip induced relaxation 
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;i=0,G=5 




FIG. 4: Fano factor versus polarization of contacts P for dif- 
ferent values of 9 in the limit of small spin-flip scattering 
intensity A — L/£sf -^ 0. 



small and the tunnel contacts have an appreciable po- 
larization. For antiparallel orientation the Fano factor is 
found to increases well above the unpolarized value due 
to a large spin accumulation in the N metal. In contrast, 
for the parallel orientation the shot noise is almost iden- 
tical with its normal state (unpolarized) value because of 
the spin accumulation suppression. At the intermediate 
angles we have found that the interplay between the spin 
accumulation precession and the suppression with mag- 
netizations angle causes that the Fano factor develops a 
minimum. The minimum shot noise can be substantially 
below the normal value depending on the relative conduc- 
tances of the N metal and tunnel contacts G/Gn and the 
polarization P. We further have shown that the spin-flip 
scattering diminishes the amplitude of the nonmonotonic 
behavior as well as the polarization effects. 

In spite of a few earlyS5iS& and recenljSl experimental 
studies devoted to spin-polarized shot noise in mag- 
netic tunnel junctions between ferrromagnets, to our 
knowledge, there have not been reports on shot noise 
measurements in metallic spin-valve systems. We expect 
that such experiments will be performed in near future, 
which will reveal the effects predicted in the present 
paper. 



in the normal metal. Applying the developed theory to 
a two terminal FNF structure we have found that the 
Fano factor has a nonmonotonic dependence on the mag- 
netizations angle, provided that the spin-flip intensity is 
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